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CO . ABSTRACT 
CO 

We examine the implications for the distribution of extrasolar planets based 
[ — | on the null results from two of the largest direct imaging surveys published to 



date. Combining the measured contrast cu rves from 22 of the star s observed with 



O 

. ^ ; the VLT NACO adaptive optics system by iMasciadri et al.l (j2005f ) , a nd 48 of the 

stars observed with the VLT NACO SDI and MMT SDI devices by iBiller et al. 
( 120071 ) (for a total of 60 unique stars: the median star for our survey is a 30 Myr 
K2 star at 25 pc), we consider what distributions of planet masses and semi- 
major axes can be ruled out by these data, based on Monte Carlo simulations of 
planet populations. We can set the following upper limit with 95% confidence: 
the fraction of stars with planets with semi-major axis between 20 and 100 AU, 
and mass above 4 Mj up , is 20% or less. Also, with a distribution of planet 
mass of 4tt oc M~ 116 in the range of 0.5-13 M Jup , we can rule out a power- 
law distribution for semi-major axis oc a a ) with index and upper cut-off 
of 18 AU, and index -0.5 with an upper cut-off of 48 AU. For the distribution 
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suggested by lCumming et al.l (120071 ). a power- law of index -0.61, we can place an 
upper limit of 75 AU on the semi-major axis distribution. At the 68% confidence 
level, these upper limits state that fewer than 8% of stars have a planet of mass 
>4 Mj up between 20 and 100 AU, and a power-law distribution for semi-major 
axis with index 0, -0.5, and -0.61 cannot have giant planets beyond 12, 23, and 29 
AU, respectively. In general, we find that even null results from direct imaging 
surveys are very powerful in constraining the distributions of giant planets (0.5- 
13 Mj up ) at large separations, but more work needs to be done to close the gap 
between planets that can be detected by direct imaging, and those to which the 
radial velocity method is sensitive. 

Subject headings: stars: planetary systems, instrumentation: adaptive optics 



Introduction 



There are currently well over 200 k nown extraso l ar pl anets, the bulk of which were 



discovered by radial velocity surveys (e.g. iButler et al.l (120061 )). While this field has initially 



been dominated by the study of the relatively easy-to-find Hot Jupiters (planets with orbital 
periods of order days), over the past several years ther e has been an incr e asing amount of 
data describing planets in larger orbits. In particular, iFischer fc Valentil (120051 ) compared 
radial velocity target stars with known planets to stars that had been monitored but did not 
show signs of planets; they concluded that about 5% of stars had planets of mass greater 
than 1.6Mj up , in orbits shorter than 4 years (within 2.5 AU). Addi tionally, they determ ined 
that planet fraction increased with the host star's metal abundance. IButler et al.l (120061 ) have 
also considered the distributions of semi-major axis and planet mass of known radial velocity 
planets, and found that both distributions are well-fit by power laws. ICumming et al.l (120071 ) 
have examined the biases of the radial velocit y technique, and found that the semi-major 
axis distribution found by iButler et al.l (120061 ). x P 1 , should be modified in light of 
the decreasing sensitivity of the radial velocity method with orbital distance, and suggest a 
power law index of -0.74 for period, instead (for solar-like stars, this corresponds to a power 
law distribution for semi-major axis where ^ oc a -0 61 ). 



A careful consideration of sensitivity of microlensing observations to planets by lGould et al 
(120061 ) suggests that for certain lensing geometries, at projected separations of ~l-4 AU, the 



lower limit for the frequency of Neptune-mass planets is 16%, making low-mass planets 
more common than giant planets in t h e inn er sol ar system (though we no te that the range 
of separations probed by iGould et al.l (120061 ) and IFischer &: Valentil (120051 ) do not precisely 
overlap, and the target star samples are not uniform between the two surveys). Additionally, 
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Gaudi et al.l (120021 ) found that from existing microlensing data, a third or less of M dwarfs 
in the galactic bulge have 1 Mj up planets in orbits between 1.5 and 4 AU, and <45% of M 
dwarfs have planets between 1 and 7 AU of mass 3 Mj up . 

One outsta nding quest i on is how the abundance of planets varies as one considers planets 
in longer orbits. iRaymondl (120061 ) has studied the dynamics of terrestrial planet formation in 
systems with giant planets, and found from numerical simulations that giant planets impede 
the formation of earth-like planets when the giant planet orbits within 2.5 AU, and that water 
delivery to a terrestrial planet is only possible in significant amounts when the giant planet 
is beyond 3.5 AU. The full extent to which giant planets impede (or encourage) water-rich 
terrestrial planet formation is still unknown. A greater understanding of the distribution of 
giant planets is a precursor to investigating the conditions under which habitable terrestrial 
planets form and evolve. 

The global distr i butio n of giant planets has also been considered from the theoretical 
direction. Ilda fc Linl (120041 ) have produced distributions of planets forming in disks by core 
accretion, showing a continuation of a power law from the radial velocity regime (within 2.5 
AU) for giant planets, out to about 10 AU, then trailing off at larger radii. It is possible 
that the lack of outer planets in these simulations may be due (at least in part) to the fact 
that these models do not consider the effects of planet-planet scattering after planets are 
formed, or it may simply be a function of the initial conditions of the simulation. In order to 
constrain such models it is necessary to measure the distribution of giant planets in longer 
orbits, so as to fully sample parameter space. 

With the advent of adaptive optics (AO) systems on large (~8m) telescopes, the abil- 
ity to detect and characterize planets by directly imaging the companion is becoming in- 
creasingly viable. Already planetary mass companions (in most cases ~13 Mj up at 40-300 
AU, or even lower mass objects with brown dwarf hosts) have been detected in certain 



favorable circumstan ces (e.g. co mpanions to 2MASS12 7: IChauvin et a 
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Itoh et al.l (120051 )). and numerous 



surveys are underway for planets around nearby, young stars (since a self-luminous planet 
is brightest at young ages). While the paucity of traditional planets (that is, planets <13 
Mj up and <40 AU orbiting a star) detected by this method has been disappointing, in this 
paper we consider how even a null result from these direct imaging surveys can be used to set 
constraints on the population of giant planets. As the sensitivity of radial velocity surveys 
to planets at larger separations decreases (due both to the smaller radial velocity signal, and 
the much longer orbital period requiring a longer time baseline of observations to adequately 
constrain the orbital parameters), at orbits wider than 10 AU only direct imaging is efficient 
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at characterizing the extrasolar planet population. 



Janson et al.l (120071 ) used VLT SDI data (part of the data considered in this work) 



for the known planet host star e Eridani, to se arch for the radial ve locity planet, given its 
predicted position from the astrometric orbit of Benedict et al.l (120061 ). Though upper limits 
were found for the planet of Mh ~19, the predicted flux of the planet could be up to 10 
magnitudes fainter, given the likely age of the system of 800 Myr. While this is young 
compared to the rest of the radial velocity planet host stars, it is quite old by the standards 
of direct imaging planet searche s, so the inabil i ty to detect this planet's flux is unsurprising. 
Previous work has been done by lKasper et al.l (120071 ) to study the region of parameter space 
unprobed by the radial velocity method, large orbital separations, by observing 22 young, 
nearby stars in the L-band from the VLT. The null result from this survey was used to set 
constraints on combinations of power-law index and upper cut-off for the distribution of the 
observed separation (not semi-major axis) of extrasolar planets. 



Masciadri et al.l (120051 ) conducted a survey of 28 young, nearby stars, with a null result 



for planets. They found that their observations were sufficient to detect a 5 Mj up planet 
at projected separations greater than 14 AU around 14 of their target stars, and above 65 
AU for all 28 stars. Similarly, their observations would have been sensitive to a 10 Mj up 
planet with a projected separation of 8.5 AU or beyond for half their sample, and greater 
than 36 AU for the full sample. These results (obtained b y adopting pu b lished ages for the 
target stars, and using the appropriate planet models of iBaraffe et al.l (120031 )) point to a 
rarity of giant planets at large separations from their parent stars. In this paper we enlarge 
the sample, and consider the implications of our null result with respect to the full orbital 
parameters of potential planets. We aim to set quantitative limits on the distribution of 
planets in semi-major axis space, and statistically rule out models of planet populations. 



2. Observations 



We begin with contrast plots (sensitivity to faint companions as a function of angular 
separation from the target star) from two s urvey s for extrasolar planets, using large tele- 
scopes and adaptive optics. iMasciadri et al.l (120051 ) carried out a survey of 28 young, nearby, 
late-type stars with the NACO adaptive optics system at the 8.2 meter Very Large Tele- 
scope (VLT). These observations have exposure times of order 30 minutes, with stars being 
observed in the H or Ks bands. Subsequent to these observations, a survey of 54 young, 
nearby stars of a variety of spectral type s (between A and M) was conducted between 2003 
and 2005, with the results reported in iBiller et al.l (120071 ). This second survey used the 
Simultaneous Differential Imager (SDI) at the 6.5 meter MMT and the 8 meter VLT, an 
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adaptive optics observational mode that allows higher contrasts by imaging simultaneously 
in narrow wavelength regions surrounding the 1.6 nm methane feature seen in c ool brown 
dwarfs and expected in extrasolar planets (ILenzen et al.ll2004l ; Oose et al.ll2005l ). This al- 
lows the light from a hypothetical companion planet to be more easily distinguishable from 
the speckle noise floor (uncorrected starlight), as the two will have very different spectral 
signatures in this r egion. This translates to higher sensitivity at smaller separations than 
the observations of IMasciadri et al.l (120051) . which were co nducted before the VLT SDI de- 
vice was commissioned (see Fig. 14 of iBiller et al.l (120071 ) for a more detailed comparison 
of the two surveys). For most of these SDI targets, the star was observed for a total of 40 
minutes of integration time, which includes a 33 degree roll in the telescope's rotation angle, 
in order to separate super speckles-which are created within the instrument, and s o will not 
rotate-from a physical companion, which will rotate on the sky (IBiller et al.ll2006l ). 

For both sets of target stars, contrast curves have been produced which give the 5oQ 
noise in the final images as a function of radius from the target stars, and thus an upper limit 
on the flux of an unseen planet in the given filter of the observations. As no planets were 
detected in either survey at the 5a level, we use these contrast curves to set upper limits on 
the population of extrasolar planets around young, nearby stars. 



2.1. Target Stars 



We construct a tar get list using 22 sta rs from the IMasciadri et al.l (120051 ) survey, and 48 
stars from the survey of IBiller et al.l (120071 ). for a total of 60 targets (10 stars were observed 
by both surveys). This first cut was made by considering stars from the two surveys that 
had contrast curves, and stars whose age could be determined by at least one of: group 
membership, lithium abundance, and the activity indicator R'hk (in three cases, ages from 
the literature were used, though these are stars that are generally older than our sample as 
a whole, and so uncertainties in the assumed ages will not adversely affect our results). Ages 
are determined by taking the age of the moving group to which the target star belongs; if the 
star does not belong to a group, the lithium or R'hk age is used, or the two are averaged if 
both are available. Lithium ages are found by com paring to lithium abundances of members 
of clusters of known ages, and similarly for R'hk (IMamajekl 120071 ). We give the full target 
list in Tabled! and details on the age determination in Table [2J We also plot our target stars 
in Fig. [TJ Overall, our median survey object is a 30 Myr K2 star at 25 pc. 



1 We note that for the SDI o bservations this th reshold corresponds to independent 5<r measurements in 
both the 0° and 33° images, see IBiller et al.l ( 2007h for details. 
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3. Monte Carlo Simulations 

In order to place constraints on the properties of planets from our null results, we run 
a series of Monte Carlo simulations of an ensemble of extrasolar planets around each target 
star. Each simulated planet is given full orbital parameters, an instantaneous orbital phase, 
and a mass, then the planet's magnitude in the observational band is determined from 
these properties (using the target star's age and distance, and theoretical mass-luminosity 
relations) as is its projected separation from the star. Finally, this magnitude is compared 
to the measured contrast curve to see if such a planet could be detected. Determining which 
simulated planets were detected, and which were not, allows us to interpret the null result 
in terms of what models of extrasolar planet populations are excluded by our survey results. 



3.1. Completeness Plots 



As in lBiller et al.l (120071 ). we use completeness plots to illustrate the sensitivity to planets 
as a function of planet mass and semi-major axis. To do this, for each target star, we create 
a grid of semi-major axis and planet mass. At each grid location we simulate 10 4 planets, 
and then compute what fraction could be detected with the contrast curve for that star. 

In general, most orbital parameters are given by well-known distributions. Inclination 
angle has a constant distribution in sin(i), while the longitude of the ascending node and the 
mean anomaly are given by uniform distributions between and 2n. Since contrast plots 
are given in terms of radius alone, it is not necessary to consider the argument of periastron 
in the simulations. 

To simulate the eccentricities of the planet orbits, we examine the orbital parameters 
of known extrasolar planets from radial velocity surveys. We consider the orbits of planets 



given by iButler et al.l (120061 ). and show their distribution of eccentricities in Fig. [2J By 
dividing the sample into two populations, based on a cut at an orbital period of 21 days, 
we can separate out the population of Hot Jupiters, which we expect to have experienced 
orbital circularization as a result of their proximity to their host stars. For both sets of 
populations, we fit a simple straight line to the distributions (the logarithmic bins for the 
Hot Jupiter population means this line translates to a quadratic fit). We note that the Hot 
Jupiter fit is plagued by small number statistics, and so the fit is likely to be less reliable 
than that for long period planets. Even for our closest target stars, such an orbital period 
gives star-planet separations less than 0.1", a regime where our contrast curves show we are 
not sensitive to planets. As a result, the manner in which the orbits of Hot Jupiters are 
simulated has effectively no impact on our final results. 
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For each simulated planet, the on-sky separation is determined at the given orbital 
phase, and mass is converted into abso lute H or Ks magnitude, following the mass-luminosity 
relations of both lBurrows et al.l (120031 ) and lBaraffe et al.l (120031 ); both of these sets of models 



have s how n success at predict ing the properties of young br own dwarfs (e.g. IStassun et al. 



(l2007h and ldose etaD (j20Q7ah ). In the case of the models of burrows et al.l J200J 



we use a 

Vega spectrum to convert the var i ous m odel spectra into absolute H and Ks magnitudes. We 
also note that the iBurrows et al.l (120031 ) models only cover a range of planet masses greater 
than 1 M Jup , and ages above 100 Myr. Since the range of ages of our target stars extends 
down to 2 Myr, and we wish to consider planets down to masses of 0.5 Mj up , we perform a 
simple extrapolation of the magnitudes to these lower ages and masses. While this solution 
is clearly not ideal, and will not reflect the complicated physical changes in these objects as 
a function of mass an d age, we feel that this method provides a good estimation of how the 
Burrows et al.l (120031 ) models apply to our survey. 



At this point, we use the distance to the target star, as well as its 2MASS flux in 
either H or Ks, to find the delta- magnitude of each simulated planet. With this, and the 
projected separation in the plane of the sky, we can compare each simulated planet to the 
5a contrast curve, and determine which planets can be detected, and which cannot. We 
also apply a minimum flux limit, based on the exposure time of the observation, as to what 
apparent magnitude for a planet is required for it to be detected, regardless of its distance 
from the parent star. For the SDI observations, which make use of optimized (compared to 
basic H-band observing) methane filters, we add an additional factor of AH=0.6 magnitudes 



(appro priate for a T6 spectral type, a conservative estimate for young planets; see lBiller et al. 



(120071 ) for details on this factor). Also, for these SDI observations, we place an upper cut-off 
on masses where, for the age of the system, the planet reaches an effective temperature of 
1400 K. Above this temperature, methane in the atmosphere of the planet is destroyed, and 
the methane feature disappears, so that the SDI subtraction now attenuates any planets, 
as well as stellar speckles. While non- methane objec t s furt her out than 0.2" are not totally 
removed in the image (e.g. Fig. 4 of Nielsen et al. ( 20051 )). for consistency we ignore this 
possibility when considering upper limits. 

In many cases in our survey, a single target star was observed at several epochs, in 
some cases with different observational parameters (such as VLT NACO SDI and VLT 
NACO Ks broadband) or even different telescopes (MMT and VLT). As a result, to be 
considered a null detection, a simulated planet must lie below the 5cr detection threshold at 
each observational epoch, and this threshold must reflect the appropriate contrast curve for 
the given observation. To account for this, for target stars with multiple observations, an 
ensemble of simulated planets is created for the earliest observational epoch, as described 
above, and compared to the contrast curve for that observation. The simulated planets then 
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retain all the same orbital parameters, except for orbital phase which is advanced forward 
by the elapsed time to the next observational epoch, and the simulated planets are now 
compared to the contrast curve from the later epoch (and so on for every available contrast 
curve). A planet that lies above the contrast curve at any epoch is considered detectable. 
Typically this elapsed time is about a year, and so is a minor effect for planets with long- 
period orbits; we nevertheless include this complexity for c ompleteness. The major benefit 
of t his method is that fo r stars observed both with SDI (jBiller et al.l 120071 ) and at H or 
Ks (IMasciadri et al.ll2005l ). it is possible to leverage both the higher contrasts at smaller 
separations with SDI and the insensitivity to the methane feature of broadband imaging, 
which allows planets of higher masses to be accessed. The e p ochs used when considering 
each obsery a tion a re those given in Table 3 of IMasciadri et al.l (120051 ) and Tables 2 and 3 of 
Biller et all (bo07h . 



We plot an example of this simulation at a single grid point in mass and semi-major axis 
in Fig. [31 for the target star GJ 182, the 18th best target star in our survey, using the planet 
models of iBurrows et al.l (120031 ) . 10 4 simulated planets (only 100 are plotted in this figure, 
for clarity) are given a single value of mass (6.5 Mj up ) and semi-major axis (10 AU). Since 
each planet has unique orbital parameters (eccentricity, viewing angle, and orbital phase), 
the projected separation varies from planet to planet, so some are above the 5a detection 
threshold of the contrast curve (the blue dots in the figure), while others are not (the red 
dots). For this particular target star and simulated planets of mass 6.5 Mj up and semi-major 
axis 10 AU, 20% of these planets can be detected. 

To produce a complete contour plot, we consider a full grid of mass (100 points, between 
0.5 and 17 Mj up ) and semi-major axis (200 points, between 1 and 4000 AU), running a 
simulation as in Fig. [3] at each of the 20,000 grid points. We then plot contours showing 
what fraction of planets we can detect that have a given mass and semi-major axis, in Fig. HJ 
again for the target star GJ 182. The hard upper limit is set by the methane cut-off, where 
the planet mass becomes high enough (for the age of the given target star) for the effective 
temperature to exceed 1400 K, at which point the methane feature is much less prominent in 
the planet's spectrum. Although there exists a Ks dataset for the star GJ 182, and additional 
observational epochs with SDI, for clarity we only use a single SDI contrast curve to produce 
this figure; the full dataset is used for subsequent analysis. If GJ 182 had a planet with 
mass and semi-major axis such that it would fall within the innermost contour of Fig. HI we 
would have an 80% chance of detecting it. Obviously, these plots make no statements about 
whether these stars have planets of the given parameters, but instead simply express our 
chances of detecting such a planet if it did exist. 
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3.2. Detection Probabilities Given an Assumed Distribution of Mass and 

Semi-major Axis of Extrasolar Planets 

With the large number of currently-known extrasolar planets, it is possible to assume 
simple power-law representations of the distributions of mass and semi-major axis of giant 



planet s, which allows for a more quantitative interpretation of our null result. IButler et a 



( 120061 ) suggest a power law of the form -n-j oc M 116 for mass, while ICumming et al.l (120071 ) 



use the power law ^ oc a -0 ' 61 for semi-major axis, in order to describe the distributions 
of known extrasolar planets. We make histograms for mass (Fig. [5]) and semi-major axis 
(Fig. E]) from the parameters of all currently-known extrasolar planets (parameters taken 
from the Catalog of Nearby Exoplanets, |http:/ /exoplanets.org| in May 2007). In both cases, 
the power laws do a reasonable job of fitting the data, above 1.6 Mj up and within 3.5 AU. 
For smaller planets, or longer periods, we would expect the observational biases of the radial 
velocity method to make the sample incomplete, thus accounting for th e drop-off of planets 



P 
•( 



from what would be predicted by the power law. We echo the caution of IButler et al.l ((2006) 
that these planets are drawn from many inhomogeneous samples, but we believe with the 
relatively large numbers the derived distributions are not far off from the actual distributions. 

In general, then, if one assumes that these power laws are universal to all stars, and 
that the semi-major axis power law continues to larger separations with the same index, the 
only outstanding question is to what outer limit (or "cut-off") this distribution continues 
before it is truncated. This cut-off is a term that can be uniquely well constrained by the 
null results from our survey. We return to this issue, after considering the results from our 
survey, in Section H~3l 

For the Monte Carlo simulations using these assumptions, in addition to the other orbital 
parameters, we obtain mass and semi-major axis through random variables that follow the 
given power law distributions, and again find what fraction of planets can be detected given 
the contrast curve for that particular target star. An example of this simulation, again for 
GJ 182, is given in Fig. [TJ showing that with an assumed upper limit for semi-major axis 
of 70 AU, and a power law with index -0.61, and mass power law index of -1.16 between 
0.5 and 13 Mj up , we would be able to detect 10% of t he simulated p l anets. Again, for this 



figure, we simply show the results using the models of iBurrows et al.l (120031 ) 



4. Analysis 



Having developed the tools to produce completeness plots, as well as compute the frac- 
tion of detected planets for various assumed models of semi-major axis, we proceed to com- 
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bine the results over all our target stars in order to place constraints on the populations of 
extrasolar planets from these two surveys. 



4.1. Planet Fraction 

A simplistic description of the number of planets expected to be detected is given by 
the expression 

N oba =eo 



N(a,M)= f P (a,M)P t (a,M) 



i=l 



That is, the number of planets one expects to detect at a certain semi-major axis and mass 
is given by the product of the detection probability (Pi) for a planet of that mass (M) and 
semi-major axis (a), and the fraction of stars (f p ) that contain such a planet (or "planet 
fraction"), summed over all target stars. In this treatment, we ignore two major effects: we 
assume that there is no change in the mass or separation distribution of planets, or their 
overall frequency, as a function of spectral type of the primary; we also do not consider 
any metallicit y dependence on the planet fraction. While t hese assumptions are clearly 
incorrect (e.g. Johnson et al. ( 2007 ). Fischer fc Valenti ( 20051 )). it is a good starting point 



for considering what constraints can be placed on the population of extrasolar planets. Also, 
we note that our sample includes 24 binaries, which may inhibit planet formation, though 
most of these binaries have separations greater than 200 AU. This leaves only ten binaries 
with separations in the range of likely planet orbits that might potentially contaminate our 
results. For simplicity, we leave these binaries in our sample, and we will return to this issue 
in Section H~3l 

Using the contrast curves from each of our 60 targets stars (as in Fig. 0J, we simply 
sum the fraction of detectable planets at each grid points for all of our stars. This gives the 
predicted number of detectable planets at each combination of mass and semi-major axis, 
assuming each target star has one planet of that mass and semi-major axis (f p (a, M) = l). 

More instructively, if we assume a uniform value of the planet fraction for all target 
stars, we can solve for f p . Then by assuming a particular value for the predicted number of 
planets (SPj), our null result allows us to place an upper limit on the planet fraction at a 
corresponding confidence level, since our survey measured a value of N(a,M)=0. In a Poisson 
distribution, the probability of obtaining a certain value is given by P = e~ M ^-, which for 
the case of a null result, v = 0, becomes P = e" 11 , so a 95% confidence level requires an 
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expectation value, fx, of 3 planets. We can thus rewrite Eq. [[J using SPj = 3, as 

Ua, M) < — = — (2) 

Put another way, if we expected, from our 5a contrast curves, to detect 12 planets (SP=12, 
for f p — 1), in order to have actually detected planets from our entire survey (N=0), the 
planet fraction must be less than ^ = 25% (f p <0.25), at the 95% confidence level. Doing 
this at each point in the grid of our completeness plots allows for an upper limit on the 
planet fraction as a function of mass and semi-major axis. 



We plot the contours of this upper limit in Fig. [8], using the planet models of lBurrows et al 



( 120031 ). A general result from these data is that, again at the 95% confdience level, we would 
expect fewer than 20% of stars to have planets of mass greater than 4Mj up with semi-major 
axis between 20 and 100 AU. There appears to be no "oasis" of giant planets (more massive 
than Jupiter) in long-period orbits: at the 85% confidence level, this upper limit on the 
fraction of stars with giant planets drops to less than 10%. 



We present the same plot, this time using the COND models of iBaraffe et al] (120031 ). in 
Fig. [91 As the two sets of models predict quite similar planet NIR magnitudes, the plots are 
virtually the same. The main difference between these models is that, given t he age distri 



bution of our target stars, higher mass planets appear slightly brighter in the IBaraffe et al. 



(120031 ) models, with th e trend reversing and l ower mass planets b ecoming fainter, as com- 
pared to the models of lBurrows et al.l (120031 ). iMarley et al.l (120071 ) have recently produced 
a third set of models, which globally predict lower luminosities for giant planets. Since 
synthetic spectra for these models are not currently available, we do not examine the conse- 
quences of these models here, though we discuss possible effects in Section [5j But we note 
th at while at 30 Myr a nd at 4 Mj up there is only a ~3X de crease in the luminosity predicted 
by lMarley et al.l (120071 ) compared to lBurrows et al.l (120031 ). the temperature of these objects 
is lower, therefore increasing the number of planets with methane that c an be detected using 
SDL As a result, even with the future use of the IMarley et al.l (120071 ) models, our results 
will not change dramatically, with respect to the total number of planets to which we are 
sensitive. 



4.2. Host Star Spectral Type Effects 

From the perspective of direct imaging searches for extrasolar planets, M-stars are es- 
pecially appealing: their lower intrinsic luminosity means a given achievable contrast ratio 
allows fainter companions to be detected, and so makes the detection of planet-mass com- 
panions seem more likely. Nevertheless, there appears to be mounting evidence that even 
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if the fraction of stars with plan e ts do es not decline when moving to later spectral types 
(and the work of IJohnson et al.l (120071 ) suggest s this fraction does decrease for M stars ), 



the mean planet mass is likely to decrease (e.g. iButler et al.l (120041 ). iBonfils et al.l (120051 )). 
While it seems natural that the initial mass of the circumstantial disk (and so the mass of 
formed planets) should scale with the mass of the parent star, such a relation is not easily 
quantified for planets at all orbital separations. Additionally, it is problematic for us to 
model planet distributions for M star hosts on radial velocity planets, when these planets 
are almost entirely in systems with a host star of spectral type F, G, or K. 

In order to investigate this effect, we divide our stars by spectral type, then recompute 
what limits we can set on the planet fraction. In Fig. [10] and [11] we plot the upper limit 
on the planet fraction for only the solar-like stars (K or earlier) in our survey (45 of our 60 
target stars, this includes the one A star in our survey, HD 172555 A). As we would expect, 
the statistics in the inner contour remain largely the same, but the contours move upward 
and to the right, as less massive and closer-in planets become harder to detect against the 
glare of earlier-type stars. 



We also consider the fifteen M stars in our sample, in Fig. [T2]and [T3J The effect of the 
smaller number of stars is apparent, thou g h the shape of the contours is again roughly the 
same. If, as is suggested by I Johnson et al.l (120071 ). giant planets are less common around low 
mass stars, or less massive stars harbor less massive planets, it becomes difficult to probe 
the population of M star planets with surveys such as these. 



4.3. Constraining the Semi-Major Axis Distribution 



We now consider what constraints can be placed on planet populations if we assume a 
basic form to the distributions. In partic ular, if we take the mass power-law from currently- 
known extrasolar planets, oc M~ L16 ( Butler et al.ll2006l ). we can constrain what types of 
power laws for semi-major axis are allowed by our survey null result. To accomplish this, we 
simulate planets using a grid of power law indices and upper cut-offs for semi-major axis for 
each of our target stars. Then, the sum of the detection fractions over the entire survey gives 
the expected number of detected planets, assuming each star has one planet (for example, 
if for 10 stars, we had a 50% chance of detecting a planet around each star, we'd expect to 
detect 5 planets after observing all 10 stars). Since we've set the distribution of planets, we 
can determine the actual planet fraction: radial velocity surveys tell us this value is 5.5% for 
planets more massive than 1.6 M j vp and with periods shorter than 4 years (closer-in than 
2.5 AU) (IFischer fc Valentil 120051 ) . We can then use the mass and semi-major axis power 
laws to find the planet fraction for planets down to 0.5 Jupiter masses and out to the given 
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semi-major axis cut-off, while always preserving the value of 5.5% for the planet fraction 
for planets >1.6 Mj up and <2.5 AU. Then, by multiplying this planet fraction by the sum 
of detection probabilities, we find the expected number of planets we'd detect given each 
distribution. At this point, we can again use the Poisson distribution to convert this to a 



confidence level (CL) for rejecting the model, given our null result: CL = 1 
is the expected number of planets for that model. 



e M , where \i 



Since stellar multiplicity is likely to disrupt planet formation, we exclude all known 
stellar binaries from our target list with projected separations less than 200 AU. Since our 
results deal mainly with the inner 100 AU around our target stars, binaries that are any closer 
would greatl y influence the format i on of planets at these radii, creating an entirely different 
population. iBonavita &: Desideral (120071 ) have shown that while the over all planet fraction 
(for ra dial velocity planets, as taken from the volume limited sample of iFischer fc Valenti 
( 120051 )) is similar between single stars and wide binaries, it decreases for stars in tight binary 
systems. Our inner cut-off o n binary separation is at a larger separation than that noted in 
Bonavita fc Desideral (120071 ). but we consider planets in much wider orbits than those de- 



tectab le wit h the radial velocity metho d. Additionally, it has been shown by lQuintana et al. 



( 120021 ) and iHolman fc Wiegertl (119991 ) that terrestrial planets could form and survive in 
the a Cen AB syst e m, d espite the relatively tight (23 AU), high eccentricity (0.5) orbit. 
Holman fc Wiegertl (119991 ) also found that for most planet is stable in a binary 

system if its orbital radius is less than ~10-20% of the binar y separation. Applying this ad- 
ditional conditi on to our sample, we remove 1 star from the iMasciadri et al.l (120051 ) survey, 
and 9 from the iBiller et al.l (120071 ) sample, leaving 50 stars in our sample. We give further 
details on the binaries in our sample in Table [31 

In Fig. HH and [151 we plot the confidence with which we can reject the model for various 
combinations of power law index and upper cut-off for the semi-major axis distribution. For 
the favored model of a power law distribution given by ^ oc a~ 0M , we can place, at the 95% 
confidence level, an upper-lim it on the semi-major axis cut-off of 75 AU (94 AU using the 
models of iBaraffe et al.l (120031 ) instead of those of iBurrows et al.l (120031 )). In other words, 
if the power law index has a value of -0.61, there can be no planets in orbits beyond a=75 
AU at the 95% confidence level (29 AU at the 68% confidence level). In Fig. [161 we show 
how these assumptions of power law index compare with the distributions of known radial 
velocity planets, as well as to what confidence we can exclude various models. 
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4.4. Testing Core Accretion Models 



We also consi der more sophist icated models of planet populations, namely the core 
accretion models of llda fc Lin! (120041 ) . Using their Fig. 12, we extract all the non-Hot- Jupiter 
giant planets, and of the 200-300 resulting planets, we run our Monte Carlo simulation by, for 
each simulated planet, randomly selecting one planet from this figure, adopting its values of 
mass and semi-major axis, then ass igning it the othe r orbital elements as usual. We consider 
each of the three cases modeled by llda fc Linl (120041 ) . 



In Fig. [T7] we plot the predicted number of planets detected from these three distri- 
butions. Again, the pla net fraction for each curve is set to match the planet fraction of 



Fischer fc Valentil (120051 ) for planets above 1.6 Mj up and within 2.5 AU. Since the predicted 
total number of planets detected range between about 0.6 and 0.7 at the end of our survey, 
we can not place any stro ng constraints on these models from our null result. For the three 
cases of llda &: Linl (120041 ). A, B, and C, we can only "rule them out" at the confidence levels 
of 45%, 49%, and 50% respectively, and again only after leaving all binaries in the sample. 
Additionally, since we are considerin g target stars of all spectral type, we are not staying 



faithful to the original simulations of llda fc Linl (12004), which consider only solar mass host 



stars. In summary, the core-accretion simulations of 
with our results. 



Ida fc Linl (120041 ) are quite consistent 



5. Discussion: Systematic Effects of Models on Results, and Other Work 



We und erscore the dep e ndenc e of th ese results upon the accuracy of the mass-luminosity 
relations of iBurrows et al.l (120031 ) and iBaraffe et al.l (120031 ) . In particular, these models 
utilize the "Hot Start" method for giant planet formation, at odds with the core ac cretion 
mechanism suggested b y the planet-metalli city relation of iFischer fc Valentil (120051 ) . The 
giant planet models of iMarley et al.l (120071 ) incorporate formation by core accretion, and 
predict systematically fainter fluxes for these young planets (typically ~3 times fainter for a 
30 Myr, 4 Mj up planet, yet the overall effect is difficult to predict without detailed models 
and spectra). Another result of moving to these models, however, would be that these planets 
are also cooler, so that the SDI method (limited to objects with effective temperatures lower 
than 140 K) will likely reach pla nets of higher masse s than would be predicted by the 
models of IBurrows et al.l (120031 ) and lBaraffe et al.l (120031 ) . 



It is possib le to e nvisio n a scenario with extrasolar planets being built by both disk 
instability (e.g. iBosd (120071 )) and core accretion, with the two types of planets segregated 
in orbital distance: inner planets being more common in orbit around metal-rich stars, 
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consistent with core accretion, while outer planets (the type to which the surveys discussed 
here are sensitive) form by disk instability. In that case, the use of the Hot Start models would 
be entirely reasonable, as these models have been shown to be mostly c onsistent wi t 



i youn g 



low-mass objects that likely form in this way (e.g. IStassun et al.l (120071 ). IClose et al.l (l2007al )). 
This possibility (which we again note is pure speculation) endangers any conclusions drawn 
from Fig. [TH and [151 which assume a single, consistent population of planets, not allowing for 
the possibility of two overlapping populations (such as one described by broken power laws). 
Our results for the upper limit on planet fraction would remain valid, however, since these 
make no assumptions on extrasolar planet populations beyond the eccentricity distribution 
(a minor factor) and the mass-luminosity relation. 

Clearly, these constraints would be stronger with a larger sample size to improve our 
statistics. Such an increase in sample size is hampered by the limited number of young, 
nearby stars: observing older targets tends to require an order of magnitude increase in 
number of targets so as to assure a similar number of detected planets. The greatest im- 
provement in these results is likely to come with more advanced planet-finding techniques, 
which increase the contrast and inner working angle to which one can detect planets close 
to their parent stars. Two such Extreme AO systems, slated to come online in the next 
several years, are VLT-SPHERE and the Gemini Planet Imager (GPI). Surveys of a sample 
of young, nearby stars (likely very similar to the target list of this work) with these planet- 
finders should be able to greatly close the gap between the sensitivities to planets of direct 
imaging and radial velocity surveys. 

While radial velocity surveys continue to have great success in finding planets, the 
limiting factor is orbital time: a planet at 10 AU takes over 30 years to complete a single 
orbit, and radial velocity planets are generally not confirmable until at least one orbit has 
elapsed. As a result, the onus on determining the characteristics of giant planets beyond 
~10 AU is largely upon direct imaging surveys. 

A survey planned for the immediate future uses the NICI (Near Infrared Coronographic 
Imager) instrument currently being commissioned on the Gemini South Telescope, with plans 
for a 50-night survey for extrasolar giant planets. It is hoped, of course, that these future 
surveys will produce actual detections, not just more null results, which when considered 
alongside the targets that were not found to harbor planets, should continue to constrain 
parameter space on the distribution of outer extrasolar giant planets. 

Another direct imaging survey for giant planets has recently been comp leted, searching 



for co mpanions to 79 young, nearby stars: the Gemini Deep Planet Survey (ILafreniere et al. 



2007 1 . For comp l etene ss, we run an extra set of simulations to compare our results to theirs. 
Lafreniere et al.l (120071 ) consider the case of planets with masses between 0.5 and 13 Mj up , 
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governed by a power law of index -1.2 (quite similar to our value of -1.16), and with a power 
law of index -1 for semi-major axis. They then set an upper limit on the planet fraction in 
three ranges of semi-major axis: 28% for 10-25 AU, 13 % for 25-50 AU, and 9.3% for 50-200 
AU, all at the 95% confidence level, using the models of lBaraffe et al.1 (120031 ). Adopting these 
same simulation parameters, we find upper limits on planet fractions of 37%, 24%, and 28%, 
respec tively. We attribute ou r somewhat lower sensitivity to the increased number of stars 
in the lLafreniere et al.1 (120071 ) survey, as well as their increased field of view (9" compared to 
the 2.2" for SDI), which makes their method bett er-suited to dete c ting p lanets at the very 
large orbital radii of the last two bins. Also, the Lafreniere et al. (j2007T ) survey was more 



consciously focused on closer stars: all 85 of their target stars are within 35 parsecs, 18 of our 



60 sta rs are beyond 35 pc. The overall results of both our work and that of lLafreniere et al. 



( 120071 ). however, are in good agreem ent for the case o f plane ts in shorter orbits: for example, 
we reach the same upper limits as lLafreniere et al.1 (120071 ) reached at the 95% confidence 
level, if we degrade our confidence level to 89% for 10-25 AU, 80% for 25-50 AU, and 63% 
for 50-200 AU. Hence the conclusions from both papers are the same: giant planets are rare 
at large separations. 

We also note that the value of the planet fraction in these intervals can be estimated 
from the uniform detectability sample of iFischer fc Valentil (120051 ). which gives 5.5% of stars 
having planets within 2.5 AU, and more massive than 1.6 Mj up . When using a model of 
planet mass with index -1.2, and semi-major axis power law index -1, as above, the planet 
fractions for the semi-major axis bins 10-25 AU, 25-50 AU, and 50-200 AU become 2.1% , 
1.6 %, and 3.2%, resp e ctivel y. It should be noted that the samples of lFischer fc Valentil (120051 ) 
and lLafreniere et al.1 (120071 ) (a s well as the one d i scusse d in this paper, for that matter) are 
not directly comparable, as the lFischer fc Valentil (120051 ) sample is made up primarily of older 



stars (>1 Gyr), and exclusively FGK spectral types, whereas the sample of lLafreniere et al. 



(120071 ) is made up of younger stars, and contains stars of M spectral type. These two effects 
push the planet fractions in opposite directi ons: younger stars are m ore likely to be metal- 
richB, and so have a higher plane t fraction (Fischer fc Valenti 2005 ). whereas M stars are 
less likely to harbor giant planets (I Johnson et al.ll2007l ). Overall, then, the upper limits from 
both papers are consistent with the predictions from radial velocity detections, with respect 
to this particular model of planet populations. 



Finally, we note that although four of our target stars do, in fact, harbor extrasolar plan- 



2 Although Table 1 of lLafreniere et al.1 (|2007l ) gives the metallicity for most of their target stars, wh ich 
give a median value of [Fe/H] = 0, more metal poor than the overall sample of Fischer fc Valentil ( 2005 ) by 
~0.1 dcx, it is notoriously difficult to make accurate metallicity measurements of young stars. As a result, 
it is likely that these reported metallicities are systematically lower than their actual values. 
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ets (HIP 30034 (AB Pic) has a wide (5.5") companion at the planet/brown dwarf boundary, 
while Eps Eri, HD 81040, and HD 128311 all have radial velocity planets), our survey can 
be regarded as a null result. Even though these planets were orbiting our target stars, we 
were unable to detect them, as they were either outside our field of view (as with AB Pic 
B), or too faint (due to their host star's age) to be detected from our images, as was the 
case with the radial velocity planets. The motivation behind our simulations is to find what 
population of hidden (undetected) planets are consistent with a lack of planet detections, 
and the knowledge of existing planets around some target stars does not change this. 



6. Conclusion 



Even without detecting extrasolar planets from our surveys, the null results provide a 
basis for setting lim its on the allowable d istribution of giant planets. From our data, using 
the planet models of lBurrows et al.1 (120031 ). we can exclude any model for planet distributions 
where more than 20% of stars of all spectral types have planets more massive than AMj up 
between 20 and 100 AU, at 95% confidence (this upper limit becomes 8% of stars with 
such planets at the 68% confidence level). If we create simple models of planet populations 
with the semi- major axis distribution governed by the power law ^ oc a", and mass by 



dN 



CX M 



-1.16 



(M/ _-v ', we can exclude giant planets in the case of a = b e yond 18 AU, and with 

a = — 0.5 beyond 48 AU. Using the distribution of ICumming et al.1 (120071 ). based on radial 
velocity observations, with a = —0.61, there can be no giant planets beyond 75 AU. All 
these statements are at the 95% confidence level; for the 68% confidence level, these upper 
limits for the outer cut-offs of giant planets become 12 AU, 23 AU, and 29 AU, for power 
law indice s of 0, -0-5 , and -0.61, respectively. With our data, the most we can say of the 
models of llda fc Lin! (120041 ) is that they are consistent with our observations at the ~50% 
confidence level. We again note that these conclusions are highly dependent on the models of 
planet luminosity as a function of the planet's age and ma ss. Additionally, we caut ion that 
since our sample differs from the volume-limited sample of iFischer fc Valentil (120051 ). known 
correlations of planet fraction with stellar mass and metallicity will likely shift our results 
from the values reported here. Nevertheless, the analysis presented here is an important first 
step in constraining the populations of extrasolar giant planets. 
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Table 1. Target Stars 



Target 




RA 1 




Dec 1 


Distance (pc) 2 


Sp. Type 


Age (Myr) 


V 1 


H 3 


Ks 3 


Obs. Mode 4 


Biller et al. (2007) 


HIP 1481 


00 


18 


26. 


1 


-63 28 39.0 


40.95 


F8/G0V 


30 


7.46 


6.25 


6. 


15 


VLT SDI 


HD 8558 


01 


23 


21. 


2 


-57 28 50.7 


49.29 


G6V 


30 


8.54 


6.95 


6 


.85 


VLT SDI 


HD 9054 


01 


28 


08. 


7 


-52 38 19.2 


37.15 


K1V 


30 


9.35 


6.94 


6. 


.83 


VLT SDI 


HIP 9141 


01 


57 


48. 


9 


-21 54 05.0 


42.35 


G3/G5V 


30 


8.11 


6.55 


6 


.47 


VLT SDI 


BD+05 378 


02 


41 


25 


9 


+05 59 18.4 


40.54 


M0 


12 


10.20 


7.23 


7 


07 


VLT SDI 


HD 17925 


(12 


52 


32 


1 


-12 46 11.0 


10.38 


K1V 


115 


6.05 


4.23 


4 


.17 


VLT SDI 


Eps Eri 


03 


32 


55. 


8 


-09 27 29.7 


3.22 


K2V 


800 


3.73 


1.88 


1 


.78 


VLT SDI 


V577 Per A 


03 


33 


13 


5 


+46 15 26.5 


33.77 


G5IV/V 


70 


8.35 


6.46 


6 


.37 


MMT SDI 


GJ 174 


04 


41 


18 


9 


+ 20 54 05.4 


13.49 


K3V 


160 


7.98 


5.31 


5 


.15 


VLT SDI 


GJ 182 


04 


59 


34. 


8 


+01 47 00.7 


26.67 


MlVe 


12 


10.10 


6.45 


6 


.26 


VLT SDI/Ks 


HIP 23309 


OS 


00 


47 


1 


-57 15 25.5 


26.26 


M0/1 


12 


10.09 


6.43 


6 


24 


VLT SDI/Ks 


AB Dor 


OS 


28 


44 


8 


-65 26 54.9 


14.94 


K1III 


70 


6.93 


4.84 


4 


.69 


VLT SDI 


GJ 207.1 


OS 


33 


44 


8 


+01 56 43.4 


16.82 


M2.5e 


100 


9.50 


7.15 


6 


.86 


VLT SDI 


UY Pic 


OS 


36 


56. 


8 


-47 57 52.9 


23.87 


K0V 


70 


7.95 


5.93 


5 


.81 


VLT SDI 


AO Men 


06 


18 


28. 


2 


-72 02 41.4 


38.48 


K6/7 


12 


10.99 


6.98 


6 


.81 


VLT SDI/Ks 


HIP 30030 


06 


19 


08. 


1 


-03 26 20.0 


52.36 


GOV 


30 


8.00 


6.59 


6 


.55 


MMT SDI 


HIP 30034 


06 


19 


12 


9 


-58 03 16.0 


45.52 


K2V 


30 


9.10 


7.09 


6 


.98 


VLT SDI 


HD 45270 


06 


22 


30 


9 


-60 13 07.1 


23.50 


G1V 


70 


6.50 


5.16 


5 


.05 


VLT SDI 


HD 48189 A 


06 


38 


00. 


4 


-61 32 00.2 


21.67 


G1/G2V 


70 


6.15 


4.75 


4 


54 


VLT SDI 


piOl UMa 


08 


39 


11 


7 


+65 01 15.3 


14.27 


G1.5V 


210 


5.63 


4.28 


4 


.17 


MMT SDI 


HD 81040 


09 


23 


47. 


1 


+20 21 52.0 


32.56 


GOV 


2500 


7.74 


6.27 


6 


.16 


MMT SDI 


LQ Hya 


09 


32 


25 


6 


-11 11 04.7 


18.34 


K0V 


13 


7.82 


5.60 


5 


45 


MMT/ VLT SDI/Ks 


DX Leo 


09 


32 


43 


7 


+ 26 59 18.7 


17.75 


K0V 


115 


7.01 


5.24 


5 


.12 


MMT/ VLT SDI 


HD 92945 


10 


43 


28. 


3 


-29 03 51.4 


21.57 


K1V 


70 


7.76 


5.77 


5 


66 


VLT SDI 


GJ 417 


11 


12 


32. 


4 


+35 48 50.7 


21.72 


GOV 


115 


6.41 


5.02 


4 


96 


MMT SDI 


TWA 14 


11 


13 


26 


5 


-45 23 43.0 


46. 00 5 


M0 


10 


13.00 


8.73 


8 


4!) 


VLT SDI 


TWA 25 


12 


15 


30 


8 


-39 48 42.0 


44.00 s 


M0 


10 


11.40 


7.50 


7. 


.31 


VLT SDI 


RXJ1224.8-7503 


12 


24 


47. 


3 


-75 03 09.4 


24.17 


K2 


16 


10.51 


7.84 


7 


71 


VLT SDI 


HD 114613 


13 


12 


03. 


2 


-37 48 10.9 


20.48 


G3V 


4200 


4.85 


3.35 


3 


.30 


VLT SDI 


HD 128311 


14 


36 


00. 


6 


+09 44 47.5 


16.57 


K0 


630 


7.51 


5.30 


5 


.14 


MMT SDI 


EK Dra 


14 


39 


00 


2 


+64 17 30.0 


33.94 


GO 


70 


7.60 


6.01 


5 


.91 


MMT SDI 


HD 135363 


15 


07 


56. 


3 


+76 12 02.7 


29.44 


G5V 


3 


8.72 


6.33 


6 


.19 


MMT SDI 


KW Lup 


15 


45 


47 


6 


-30 20 55.7 


40.92 


K2V 


2 


9.37 


6.64 


6 


40 


VLT SDI 


HD 155555 AB 


17 


17 


25 


5 


-66 57 04.0 


30.03 


G5IV 


12 


7.20 


4.91 


4 


70 


VLT SDI/Ks 


HD 155555 C 


17 


17 


27 


7 


-66 57 00.0 


30.03 


M4.5 


12 


12.70 


7.92 


7 


63 


VLT SDI/Ks 


HD 166435 


18 


09 


21 


4 


+ 29 57 06.2 


25.24 


GO 


100 


6.85 


5.39 


5 


.32 


MMT SDI 


HD 172555 A 


18 


45 


26 


9 


-64 52 16.5 


29.23 


A5IV/V 


12 


4.80 


4.25 


4 


.30 


VLT SDI 


CD -64 1208 


18 


45 


37. 





-64 51 44.6 


34.21 


K7 


12 


10.12 


6.32 


6 


10 


VLT SDI/Ks 


HD 181321 


19 


21 


29 


8 


-34 59 00.5 


20.86 


G1/G2V 


160 


6.48 


5.05 


4 


93 


VLT SDI 


HD 186704 


1!) 


45 


57. 


3 


+04 14 54.6 


30.26 


GO 


200 


7.03 


5.62 


5 


52 


MMT SDI 


GJ 799B 


20 


41 


51 


1 


-32 26 09.0 


10.22 


M4.5e 


12 


11.00 


0.00 


-99.00 


VLT SDI/Ks 


GJ 799A 


20 


41 


51 


2 


-32 26 06.6 


10.22 


M4.5e 


12 


10.25 


5.20 


4.94 


VLT SDI/Ks 


GJ 803 


20 


45 


09. 


5 


-31 20 27.1 


9.94 


MOVe 


12 


8.81 


4.83 


4.53 


VLT SDI/Ks 


HD 201091 


21 


06 


53. 


9 


+38 44 57.9 


3.48 


K5Ve 


2000 


5.21 


2.54 


2 


25 


MMT SDI 


Eps Indi A 


22 


03 


21. 


7 


-56 47 09.5 


3.63 


K5Ve 


1300 


4.69 


2.35 


2 


24 


VLT SDI 


GJ 862 


22 


29 


15 


2 


-30 01 06.4 


15.45 


K5V 


6300 


7.65 


5.28 


5 


11 


VLT SDI 


HIP 112312 A 


22 


44 


57. 


8 


-33 15 01.0 


23.61 


M4e 


12 


12.20 


7.15 


6.93 


VLT SDI 


HD 224228 


23 


56 


10 


7 


-39 03 08.4 


22.08 


K3V 


70 


8.20 


6.01 


5.91 


VLT SDI 


Masciadri ct al. (20051 


HIP 2729 


00 


34 


51 


2 


-61 54 58 


45.91 


K5V 


30 


9.56 


6.72 


(i 


53 


VLT Ks 


BD +2 1729 


06 


18 


28. 


2 


-72 02 42 


14.87 


K7 


30 


9.82 


6.09 


5 


.87 


VLT H 


TWA 6 


07 


39 


23. 





02 01 01 


77.00 5 


K7 


30 


11.62 


8.18 


8 


04 


VLT Ks 


BD +1 2447 


10 


18 


28. 


8 


-31 50 02 


7.23 


M2 


12 


9.63 


5.61 


5 


.31 


VLT H 


TWA 8A 


10 


28 


55. 


5 


00 50 28 


21.00 s 


M2 


115 


12.10 


7.66 


7 


.43 


VLT Ks 


TWA 8B 


11 


32 


41 


5 


-26 51 55 


21.00 s 


M5 


100 


15.20 


9.28 


9 


.01 


VLT Ks 


TWA 9A 


11 


32 


41 


5 


-26 51 55 


50.33 


KS 


800 


11.26 


8.03 


7 


.85 


VLT Ks 


TWA 9B 


11 


48 


24 


2 


-37 28 49 


50.33 


Ml 


70 


14.10 


9.38 


9 


.15 


VLT Ks 


SAO 252852 


11 


48 


24 


2 


-37 28 49 


16.40 s 


K5V 


160 


8.47 


5.69 


5 


.51 


VLT H 


V343 Nor 


14 


42 


28. 


1 


-64 58 43 


39.76 


K0V 


12 


8.14 


5.99 


5 


.85 


VLT Ks 


PZ Tel 


IS 


38 


57. 


6 


-57 42 27 


49.65 


KOVp 


12 


8.42 


6.49 


6 


.37 


VLT Ks 


BD-17 6128 


18 


53 


05 


9 


-50 10 50 


47.70 


K7 


70 


10.60 


7.25 


7 


04 


VLT Ks 



from the CDS Simbad service 

derived from the Hipparcos survey |Perrvman et al ] i 19971 ) 
3 from the 2MASS Survev lCutri et al.l 120031 ') 
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4 In cases were target stars were observed by both Masciadri ct al. (2005) and Billcr ct al. (2007). the star is listed in the Billcr ct al. (2007) section, 
with Obs. Mode given as "VLT SDI/Ks," for example. 

5 Distance from lSong et al.| j2003|) 

^Distance from Zuckcrman ct al ] j2001al) 
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Table 2. Age Determination for Target Stars 



Target 


Sp. Type* 


Li EW (mas)* 


Li Age (Myr) 




R-'ffif Age 


Group Membership 1 


Group Aj 


re 1 Adopted Age 


Biller et al. (20071 


HIP 1481 




129 3 


100 


-4.360* 


200 


Tuc/Hor 


30 


30 


HD 8558 


G6V 2 


205 5 


13 






Tuc/Hor 


30 


30 


HD 9054 


K1V 2 


170 s 


160 


-4.236 6 


<100 


Tuc/Hor 


30 


30 


HIP 9141 


G3/G5V 7 


181 8 


13 






Tuc/Hor 


30 


30 


BD+05 378 


M0 9 


15 10 








Pic 


12 


12 


HD 17925 


K1V 7 


194 8 


50 


-4.357 6 


200 


Hcr/Lyr 


115 


115 


Eps Eri 


K2V 11 






-4.598 6 


1300 






800 12 


V577 Per A 


G5IV/V 13 


219 13 


3 






AB Dor 


70 


70 


GJ 174 


K3V 14 


118 s 


160 


-4.066 6 


<100 






100 


GJ 182 


MlVe 10 


280 15 


12 










12 


AB Dor 


K1III 2 


267 s 


10 


-3.880 6 


<100 


AB Dor 


70 


70 


GJ 207.1 


M2.5e 16 














100 17 


HIP 23309 


MO/l 18 


294 ls 


12 


-3.893 6 




Pic 


12 


12 


UY Pic 


KOV 19 


263 s 


10 


-4.234 6 


<100 


AB Dor 


70 


70 


AO Men 


K6/7 18 


357 ls 


6 


-3.755 6 




Pic 


12 


12 


HD 45270 


G1V 2 


149 5 








AB Dor 


70 


70 


HD 48189 A 


G1/G2V 2 


145 s 


25 


-4.268 6 


100 


AB Dor 


70 


70 


HIP 30030 


GOV 20 


219 8 


2 






Tuc/Hor 


30 


30 


HIP 30034 


K2V 2 










Tuc/Hor 


30 


30 


piOl UMa 


G1.5V 21 


135 s 


100 


-4.400 22 


320 






210 


DX Leo 


KOV 21 


180 s 


100 


-4.234 6 


<100 


Her /Lyr 


115 


115 


HD 81040 


GOV 21 


24 23 


2500 










2500 


LQ Hya 


KOV 21 


247 s 


13 










13 


HD 92945 


K1V 21 


138 s 


160 


-4.393 6 


320 


AB Dor 


70 


70 


GJ 417 


GOV 24 


7g 25 


250 


-4.368 26 


250 


Hcr/Lyr 


115 


115 


TWA 14 


M0 27 


600 27 


8 






TW Hya 


10 


10 


RXJ1224.8-7503 


K2 2s 


250 2S 


16 










16 


TWA 25 


MO 9 


494 29 


10 






TW Hya 


10 


10 


HD 114613 


G3V 30 


100 31 


400 


-5. 118 6 


7900 






4200 


EK Dra 


GO 32 


212 s 


2 


-4.180 22 


<100 


AB Dor 


70 


70 


HD 128311 


KO 21 






-4.489 26 


630 






630 


HD 135363 


G5V 21 


220 s 


3 










3 


KW Lup 


K2V 30 


430 33 


2 










2 


HD 155555 AB 


G5IV 18 


205 s 


6 


-3.965 6 


<100 


Pic 


12 


12 


HD 155555 C 


M4.5 18 










Pic 


12 


12 


CD -64 1208 


K7 18 


580 ls 


5 






Pic 


12 


12 


HD 166435 


GO 34 






-4.270 22 


100 






100 


HD 172555 A 


A5IV/V 2 










Pic 


12 


12 


HD 181321 


G1/G2V 30 


131 s 


79 


-4.372 6 


250 






160 


HD 186704 


GO 35 






-4.350 22 


200 






200 


GJ 799A 


M4.5e 16 










Pic 


12 


12 


GJ 799B 


M4.5e 16 










Pic 


12 


12 


GJ 803 


MOVe 16 


51 8 


30 






Pic 


12 


12 


HD 201091 


K5Vc 16 






-4.704 6 


2000+ 
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~*~In general, we have only determined Ca h K a g es f° r stars with spectral types Kl or earlier, but in the case of these two K5 stars, we have only 
the TI'hk measur ement on wh i ch to rely for age determination. The calibration of Mt. Wilson S-indcx to H'hk f° r K5 stars (B-V ~ 1.1 mag) has not 
been well-defined ( Noycs ct al. (1984); specifically the photosphcric subtraction), and hence applying a R'hK" vs - a S c relation for K5 stars is unlikely to 
yield useful ages. Although we adopt specific values for the ages of these stars, it would be more accurate to state simply that these stars have ages >1 
Gyr. As a result, almost all simulated planets are too faint to detect around these stars, so the precise error in the age docs not significantly affect our 
final results. 
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Table 3. Binaries 
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Pounds et al. (19931 
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AB Dor 
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134 (Ba/Bb) 


Close et al. (2005) 


Binary M stars 


AB Dor 


0.15 (C) 
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Close et al. (2005) 


Very low-mass M Star 
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Cb^vrvinetaL (2005) 


Planet/Brown Dwarf 
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16.5 
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Biller et al. (2007) 
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Zuckerman et al. (2001a) 
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2100 
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13 
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36 


Wilson (1954) 


Target Star GJ 799B, M4.5 
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Song et al. (2002) 
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15.7 


200 


Poveda et al. (1994) 


HD 128898, Ap 


V343 Nor 


10 
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Song et al. (2003) 


M4.5 


BD-17 6128 


2 
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Neuhauser et al. (2002) 
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Fig. 1. — The 60 target stars from our two surveys (though five stars are too old to appear 
on this plot). These stars are some of the youngest, nearest stars known, spanning a range of 
spectral type. The size of the plotting symbol and the color is proportional to the absolute 
H magnitude of the star: a bigger, bluer symbol corresponds to a brighter and hotter star. 
The legend gives approximate spectral type conversions for main sequence stars, but we 
note that these stars have been plotted by their 2MASS H-band fluxes, and as a result their 
actual spectral type can vary from that shown in the legend. See Table [1] for more complete 
properties of these stars. The median target star is a 30 Myr K2 star at 25 pc. 
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Fig. 2. — The assumed distribution for the orbital eccentricities of extrasolar planets. The 
datapoints represen t the histograms for planets found to date with the radial velocity method 
(jButler et al.ll2006l ). with error bars as 1-sigma Poisson noise based on the number of planets 
per bin. Planets are divided to separate "Hot Jupiters," based on a period cut at 21 days; 
long period planets are divided into linear bins, short-period ones into logarithmic bins. In 
both cases, a simple linear fit is a good representation of the data. 
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Fig. 3. — The results of a single simulation of 10 4 planets aroun d the SDI target star GJ 182 
(for clarity, only 100 points are plotted here) ( jBiller et al.l 120071 ). Each planet has a mass of 
6.5 Mj up , and a semi-major axis of 10 AU. Due to various values of eccentricity, viewing angle, 
and orbital phase, the projected separation of each simulated planet departs from the semi- 
major axis, and the points smear across the horizontal direction, with projected separation 
running between and 14 AU. Planets that are above the contrast curve are detected (blue 
dots), while those below are not (red dots). In this case, 20% of these simulated planets were 
detected. By running this simulation over multiple grid points of mass and semi-major axis, 
we produce a full completeness plot, such as Fig. HI 
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Fig. 4. — A full completeness plot for the target star GJ 182. As a function of planet mass 
and semi-major axis (with grid points between 0.5 and 17 Mj up for mass, and semi-major 
axis between 1 and 4000 AU, though only the inner 210 AU are plotted here), the contours 
give the probability of detecting a planet with those parameters given the available contrast 
curve. At each grid point, 10 4 planets are simulated, as shown in Fig. [31 and the fraction 
that can be detected is returned. The left edge is strongly influenced by the shape of the 
contrast curve, while the right edge depends mainly on the projected field of view of the 
observation. The hard upper limit at 9 Mj up is set by the methane cut-off imposed by the 
SDI method, when the simulated planets exceed 1400 K and cease to have a strong methane 
signature in the spectrum. The fact that the contours do not precisely line up at this limit is 
simply a result of the interpolation used to plot the contours. Completeness plots for all 60 
survey stars are available online at j http : / / exoplanet . as . arizona.edu / aldose/ exoplanet . html | 
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Fig. 5. — The assumed mass distribution of extrasolar planets, plotted against the histogram 
of known planets detected by the radial velocity method . Throughou t this p aper we adopt 
a power law of the form 4tj oc M~ L16 , as suggested by iButler et al.l (120061 ). which does a 
reasonable job fitting the data. 
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Fig. 6. — The distributions that we consider for semi-major axis of extrasolar planets, again 
wi th the histogram of kn own radial velocity planets. We adopt the observed distribution 
of ICumming et al.l (120071 ) . with ^ oc a -0 - 61 , which is suggestive of the existence of wider 
planets, given that radial velocity surveys should be especially sensitive to hot Jupiters 
(producing an over-abundance at small separations) and less sensitive to long-period orbits 
(resulting in a decline in detected planets at larger separations). 
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Fig. 7. — 10 5 simulated planets around the SDI targ et star GJ 182, following the distribu- 
tions for mass f j xz oc -M -1 - 16 ) of iButler et al.l ( 20061 ) and semi-major axis oc a -0 ' 61 ) of 



Cumming et al.l (120071 ). with mass running from 0.5 to 13 Mj up , and semi-major axis cut 



off at 70 AU (since there is a range of eccentricities, separation can exceed the semi-major 
axis cut-off). Detected planets (blue dots) are those that lie above the contrast curve, above 
the minimum flux level, and below the methane cut-off. In this case, 10% of the simulated 
planets could be detected with this observation. Using the metric of completeness to planets 
with this mass and semi-major axis distribution, GJ 182 is the 18th best target star in our 
sample. 
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Burrows et al. (2003), All Stars 
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Fig. 8. — The upper limit on the fraction of stars with pla nets (f p ), as a f unctio n of mass 
and semi-major axis (see Eq. [2]), using the planet models of iBurrows et al.l (j2003l ). with the 
95% confidence level plotted as thin blue lines. We also plot in thicker red lines the 68% 
confidence level contours. Given the results of our survey, we would expect, for example, less 
than 20% (as indicated by the thin dashed blue line) of stars to have a planet of mass greater 
than 4 Mj up in an orbit 20 < a < 100 AU, and less than 50% of stars (the dot-dashed thin 
blue line) to have planets more massive than 4 Mj up with semi-major axes between 8 and 
250 AU, at the 95% confidence level. Also plotted in the solid circles are known extrasolar 
planets. There is still a gap between planets probed by direct imaging surveys (such as the 
ones described in this work), and those using the radial velocity method. 
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Baraffe et al. (2003), All Stars 
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Fig. 9. — The same as Fig. El but instead using the models of iBaraffe et al.l (120031 ) to convert 
between planet mass and NIR magnitudes. The COND models generally p redict brighter 



planet s for higher masses, but fainter planets at lower masses, compared to the lBurrows et al. 



(120031 ) models. Nevertheless, the two sets of models predict similar overall results. 
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Burrows et al. (2003), AFGK Stars 
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Fig. 10. — The 95% and 68% confidence u pper limit on p l anet f raction, limited only to stars 
of spectral type A through K, using the iBurrows et al.l (120031 ) models. Since with earlier 
spectral types the parent star is intrinsically brighter, it becomes more difficult to access 
planets of smaller masses or smaller separations. For AFGK stars we can only say, at the 
95% confidence level, that less than 20% of stars have M > 7M Jup planets at 30-70 AU, or 
a limit of 50% for planets with masses above 6 Mj up at 10-200 AU. 
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Fig. 11. — The same as Fig. [TOj but with the lBaraffe et al.l (120031 ) models used to find planet 
masses. 
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Fig. 12. — Now using only our 15 M st ars, we again p l ot the 95% and 68% confidence level 
upper limit on planet fraction, using the lBurrows et al.l (120031 ) models. While the plot follows 
the shape of Fig. El the removal of about three-quarters of the target stars reduces the upper 
limit that can be set on the planet fraction. Hence less than 50% of M stars should have 
planets wi th M > 4Mj „, „ from 10 to 80 AU, at 95% confidence. The analysis of microlensing 
results by iGaudi et al.l (120021 ) sets upper limits on the planet fraction of M dwarfs in the 
galactic bulge of 45% for 3 M Jup planets between 1 and 7 AU, and 33% for 1 M Jup planets 
between 1.5 and 4 AU. While even our 50% conto ur (at the 68% conf idence level) does not 
probe the area of parameter space considered by I G audi et al.l (120021 ) . which places upper 
limits on 1 Mj up planets between 1.5 and 4 AU around M dwarfs of <33%, and <45% for 3 
Mj up planets between 1 and 7 AU, the microlensing upper limits are unsurprising given our 
limits at somewhat larger separations for planets of the same mass. Though we note that the 
composition (especially in terms of stellar metallicity) is li kely to differ gr e atly b etween the 
two samples. Also, we again draw attention to the fact that iJohnson et al.l (120071 ) has shown 
that for M stars, giant planets at small radii are less common than around more massive 
stars. 
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Fig. 13. — As with Fig. [T2], only now with the 
planet masses. 
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Fig. 14. — The confidence level with which we can reject models of planet populations, 
assuming a power-law distribution for semi-major axes oc a a ), as a function of the 
power law index and upper cut-off (N(a)=0 for a > ac u t-off) - The expected powe r-law 
index from the radial velocity distribution (see Fig EJ) is -0.61 (jCumming et al.l 120071 ). and 
given these data we can place a 95% confidence limit on the upper cut-off of 75 AU. At 68% 
confidence, there cannot be giant planets in or bits beyond 29 AU, fo r this choice of power 
law index. For this figure, we use the models of iBurrows et ajj (120031 ) 
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Fig. 15. — The same as Fig. [HJ but using the models of iBaraffe et al.l (120031 ). The 95% 



confidence upper cut-off for semi-major axis for the oc a a model now moves to 94 AU 
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Fig. 16. — The histogram (in blue) of the distribution of known extrasolar giant planets 
found with the radial velocity method, plotted against a series of power laws considered in 
Fig. [HlandCLHl Since radial velocity observations are only complete to about 2.5 AU, a less 
steep drop-off of planets with semi-major axis is possible. We give the confidence with which 
we can rule out various combinations of power law index and upper cut-off (the percentages 
in red), for indices of -1, -0.61, -0.25, and upper cut-offs of 10 AU, 20 AU, 40 AU, and 80 AU. 
While we have insufficient statistics to place strong constraints on the power law of index -1, 
we can rule out the other two with increasing confidence as larger values of the upper limit 
are considered. For example, a power law of the form ^ oc a -0 ' 25 must cut-off at 26 AU 
(95% confidence), while the most likely power law of index -0.61 must have its cut-off at 75 
AU (also at the 95% confidence level). 
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Fig. 17. — The number of planets we would expect to detect at the end of the survey, as 
a function of the number of target stars observed, out of our total sample of 60. Stars are 
divided into bins based on binarity, and within each bin the stars are arranged so that the 
best targets are observed first. The first four models use power laws with oc a -061 , with 
the upper cut-off given. These models can be ruled out with increasing confidence with cut- 
offs beyon d 40 A U as increasingly close binaries are added to the sample. Since the three 



Ida fc Linl (120041 ) models predict less than one planet from our survey, we can only place 



very limited constraints on these models at this time, namely that cases A, B, and C are 
inconsistent with our null result at the 45%, 49%, and 50% confidence levels, respectively, if 
all binaries are included in our sample. 



